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In this note we prove an existence and uniqueness result for the solution of multidimensional 
stochastic delay differential equations with normal reflection. The equations are driven by a frac- 
tional Brownian motion with Hurst parameter H > 1/2. The stochastic integral with respect to 
the fractional Brownian motion is a pathwise Riemann-Stieltjes integral. 
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1. Introduction 

Consider a stochastic delay differential equation with positivity constraints. More pre- 
cisely, we deal with a stochastic delay differential equation with normal reflection on M'* 
of the form 

X(t) = 7]{0)+ [ b{s,X)ds+ [ <j{s,X{s-r))dW," + Y{t), te{0,T], 



JO 



(1.1) 

Xit)=7]{t), te[~r,0]. 

Here, r denotes a strictly positive time delay, = {W^'^,j ~ l,...,m} are inde- 
pendent fractional Brownian motions with Hurst parameter H > ^ defined in a com- 
plete probability space (17,J^, P), the hereditary term b{s,X) depends on the path 
{X{u),—r <u< s}, while ri:[—r,0] — >■ R^J. is a non-negative smooth function with 
> for i = 1, . . . , d} and F is a vector- valued non-decreasing process 
which ensures that the non-negativity constraints on X are enforced. 
Set 

Z{t)=r]{Q)+ [ b{s,X)ds+ f a(s,X(s-r))dWf , te[0,T]. (1.2) 
Jo Jo 
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It is known that we have an exphcit formula for the regulator term Y in terms of Z: for 
each i = 1, . . . ,d, 

Y\t) = max {Z\s))-, t £ [CTl. 
se[o,t] 



The solution of (1.1) then satisfies 

X{t) 



[Z{t) + Y{t), te[o,T], 
\in{t), te[-r,o]. 



We call (1.1) a delay differential equation with rcfiection with hereditary drift driven 
by a fractional Brownian motion. To the best of our knowledge, this problem has not 
been considered before in the extensive literature on stochastic differential equations. 

There are many references to stochastic systems with delay (see, as a basic refer- 
ence, [9]), but the literature concerning stochastic differential equations with delay driven 
by a fractional Brownian motion is scarce. The existence and uniqueness of solutions 
[2, 7, 10, 12], existence and regularity of the density [7] and convergence when the delay 
goes to zero [3] have all been studied. 

On the other hand, there has been little work on stochastic differential equations 
with delay and non-negativity constraints. We can only refer the reader to the book by 
Kushner [5] dedicated to the study of numerical methods for this class of equations, and 
the paper of Kinnally and Williams [6], where the authors obtain sufficient conditions 
for existence and uniqueness of stationary solutions for stochastic differential equations 
with delay and non-negativity constraints driven by standard Brownian motion. 

As described in the paper of Kinnally and Williams [6] , there are some models affected 
by some types of noise where the dynamics are related to propagation delay and some 
of them are naturally non-negative quantities, for instance, applications such as rates 
and prices in Internet models, and concentrations of ions or proportions of a population 
that are infected (see the references in [6]). Therefore, it is natural to consider stochastic 
differential equations with delay and non-negativity constraints. In this paper, we initiate 
the study when the noise is not a standard Brownian motion but a fractional Brownian 
motion. 

The main novelty of this paper is the use of non-negative constraints dealing with 
fractional Brownian motion. We have used Skorokhod's mapping. 
Set 

C+(M+,R'') :={a;eC(M+,M'*):.T(0) GM:|^}. 
We now recall the Skorokhod problem. 

Definition 1.1. Given a path z G C-|-(R+,R'^), wt say that a pair {x,y) of functions in 
C+(]R+,R'^) solves the Skorokhod problem for z with reflection if: 

1. x{t) = z{t) + y{t) for allt>0 and x{t) G for each t > 0; 

2. for each i = 1, . . . ,d, y'(0) = and y* is non- decreasing; 

3. for each i = 1, . . . , d, a;*(s) dy*(s) ~ for all t > 0, so y' can increase only when 
is at zero. 
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It is known that we have an exphcit formula for y in terms of z: for each i = 1, . . . , d, 



The path z is called the reflector of x and the path y is called the regulator of x. We 
use the Skorokhod mapping for constraining a continuous real- valued function to be non- 
negative by means of reflection at the origin. We will apply it to each path of z defined 
by (1.2). Note that because we are dealing with a multidimensional case, the mapping 
will be applied to each component. 

We must also explain how to understand the stochastic integral appearing in (1.1). 
Since H > ^, the stochastic integral in (1.1) is defined using a pathwise approach. Indeed, 
if we have a stochastic process {u{t),t > 0} whose trajectories are A-Holder continuous 
with X> 1 — H, then the Riemann-Sticltjes integral u{s) dW^ exists for each trajec- 
tory (see Young [13]). Using the techniques introduced by Young [13] and the p-variation 
norm, Lyons [8] began the study of integral equations driven by functions with bounded 
p- variation, where p € [1,2). Zahle [14] introduced a generalized Stieltjes integral using 
the techniques of fractional calculus. The integral is expressed in terms of fractional 
derivative operators and coincides with the Riemann-Stieltjes integral J^^ fdg when the 
functions / and g are Holder continuous of orders A and /3, respectively, with X + (3 > 1. 
Using this Riemann-Stieltjes integral, Nualart and Rascanu [11] obtained the existence 
and uniqueness of a solution for a class of multidimensional integral equations. 

In this paper, also using the Riemann-Stieltjes integral, we will prove the existence and 
uniqueness of a solution to equation (1.1). Our results are inspired by those in Nualart 
and Rascanu [11] and Ferrante and Rovira [3]. Using some estimates presented in those 
papers, we will first prove our results for deterministic equations and will then easily 
apply them pathwise to fractional Brownian motion. 

Since our definition of the stochastic integral holds for iJ > i , we cannot extend our 
approach to the case G (^, ^). However, in a forthcoming paper, we will use the method 
used by Hu and Nualart in [4] to consider the case H < ^. 

The structure of the paper is as follows. In the next section, we give our hypothesis 
and we state the main results of our paper. In Section 3 we give some useful estimates for 
Lcbcsgue and Riemann-Stieltjes integrals inspired by the results in [11] and [3]. Section 4 
is devoted to proving our main result: the existence, uniqueness and boundedness of so- 
lutions to deterministic equations. In Section 5 we recall how to apply the deterministic 
results to the stochastic case, while the Appendix is devoted to giving some technical re- 
sults such as a fixed point theorem and some properties related to the Skorokhod problem. 



2. Main results 

Let a G (0, ^) and r > 0. Let {s,t) C [— r, T] and denote by W^'°°(s,i;M'^) the space of 
measurable functions / : [s, i] M'' such that 



y\t)= max (z''(s)) . 
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For any < A < 1, denote by C^{s,t;M.'^) the space of A-H61der continuous functions / : 
[s,t]^ R'^ such that 

ll/lU(M):=ll/lloo(M)+ sup <oo, 

s<u<v<t ["^ 



where 



\oo(s,t) ■= sup \f{u) 

ue[s,t] 



Since we will use the spaces W(5''°°(— r, r;M'*) and C'^(— r, r;R'*) extensively, we will use 

the notation ||/|U,oo(r) := ll/IU,oo(-r,T), ll/IU,A(r) := ll/IU,A(-r,T), ll/IU(r) \\f\\x{-r,T) 
and ||/||oo(r) := ll/lloo(-r,T)- 

Consider the following hypothesis. 

(HI) a: [0,T] x M'' ^ M'' x R"' is a measurable function such that there exist some 
constants /3 > and Mo > such that the following properties hold: 

1. \a{t,x) - a{t,y)\ < Mo\x - y\ Va;,y £ R^Vt £ [0,r]; 

2. \ait,x) - a{s,x)\ < Ma\t - s\^^ Vx e R"*, Vt, s £ [0, T]. 

(H2) b:[0,T] x C{-r,T]R'^) -j- R'' is a measurable function such that for every t > 
and / e C{-r,T;R'^),b{tJ) depends only on {/(s);-r <s<t}. Moreover, there 
exists some bo G L''{0,t;M.'^) with p>2 and VA^ > there exists some Ln > such 
that: 

1. \b{t,x) - b{t,y)\ < LArsup_^<,<j - y{s)\ Vx,y such that ||x||oo(r) < N, 

\\y\\ooir)<NWtG[0,n 

2. |&(t,a;)| <ioSup_,<,<, \x{s)\+bn{t) Vt £ [0,r]. 

(H3) There exist some 7 G [0, 1] and Kq > such that 

\a{t,x) \ < A'o(l + VxgR'', yte[o,T]. 

Under these assumptions, wc arc able to prove that our problem admits a unique solution. 
Our main result reads as follows. 

Theorem 2.1. Assume that T] e W^''^{-r,0;R^), and that b and a satisfy hypothe- 
ses (HI) and (H2), respectively, with (3 > 1 — H . min{i, /?}. // a G (1 — -ff, ag) 
and p < ^ , then the equation (1.1) has a unique solution, 

X G i°(ri, P; Wo"'°°(-r, T; R'')), 

and for P-almost allujen, X(w, ■) G C^""(0,T;M''). 

Moreover, if a G (1 -iJ, ao V(2-7)/4) and (H3) /loZds, t/ien i;(||X||^ ^^^^ < ooVp> 1. 

Examples. Note that the following equations satisfy our hypothesis: 
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(a) (linear example) for any a, 6 € M, 

X{t)=r+ f X{s-r)ds+ f {aX{s ~r) +b)dWj^ + Y{t), te(0,T], 
Jo Jo 

X{t)=t + r, <e[-r,0]; 

(b) (non-linear example) 

X{t)= I cos{X{s))ds+ [ siii{s + X{s-r))dWf +Y{t), ^e(0,T], 
Jo Jo 

X{t)^t^, te[-r,0]. 

3. Preliminaries 

In this section, we give some useful estimates for Lebesgue and Riemann-Stieltjes inte- 
grals. These types of estimates were presented in the work of Nualart and Rascanu [11] 
and adapted to the delay case by Ferrante and Rovira [3]. Since our results are directly 
inspired by these works, we do not give the proofs, but instead direct the reader to these 
references. 

We will need to introduce a new norm in the space WQ'°°{s,t;W^): for any A > 1, 

||/IU,A(.,t) sup cxp(~A^)f|/(^.)|+ r ~ d^^ 
ue[s,t] \ Js l"" ^) J 

It is easy to check that for any A > 1, this norm is equivalent to ||/||a,oo(s,t)- 
3.1. Lebesgue integral 

Wc first consider the ordinary Lebesgue integral. Given a measurable function / : [— r, T] 
M'', we define 

Ab)i 



F^^um^ / b{uj)du. 

Jo 

We first recall some estimates that constitute an obvious adaptation of [3], Proposi- 
tion 2.2. 

Proposition 3.1. Assume that b satisfies (H2) with p = — and [s,t] C [0,r]. If f € 

Wo'°°{-r,t;R'^), then F'-''\f){-) = J^b{uJ)dueC^-"{s,t;R'^) and 

for all A > 1, where i G {1,2}, are positive constants depending only on a, t, Lq and 
Bo,a = 1|6o||li/o- 
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3.2. Riemann— Stieltjes integral 



Let us now consider the Ricmann-Stieltjes integral introduced by Zahle, which is based 
on fractional integrals and derivatives. We refer the reader to the paper of Zahle [14] and 
the references therein for a detailed presentation of this generalized Stieltjes integral and 
the associated fractional calculus. Here, we will just recall some basic facts. 

Fix a parameter < a < |. Denote by W^~"'°°(0,r;IR) the space of measurable func- 
tions 5 : [0, T] M such that 



|ff||l-a,oo,T : 



sup 



o<.<t<TV (t-s)i-" J, (y-s) 
Moreover, if g belongs to H/y~"'°°(0, T; R), the n we define 

1 



Aa(.9) 



< 



r(l -a) o<s<t'<T'^ * 

1 



sup \iDll'^9t-M\ 

<T 

3lll-a,oo,T < OO, 



r(l -a)r(a) 

where D^I°'gt- denotes a fractional Weyl derivative, 

9t- (s) = {g{s) - 5(t-))l(o,t) (s), 

where g{t~) ~ lime\,o3(^ ~ e) and F is the Euler function. We also denote by W(5''^(0, 
T; M) the space of measurable functions / on [0, T] such that 



I/Ik 



1/(^)1 



ds- 



T ps 



{s - 



duds < OO. 



Note that if / is a function in the space Wq'^{Q, T; M) and g £ l¥y~"'°°(0, T; M), then the 
integral Jq f dg exists for aU t G [0,T] and we can define 



G(/)W:= r.f(s)dg.= r/(s)l(o.t)(s) 
Jo Jo 



dgs- 
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Furthermore, the following estimate holds: 

rt 



fdg 



<Aa(.9)||/|| 



(3.1) 



Moreover, if / G Wq'°°{0,T), it is proved in [11], Proposition 4.1, that for each s <t, 



fdg 



<A„(g)C„,T(i-s)i-"l|/l|„, 



(3.2) 
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Let us consider the term 

Gl'^Hm)^ faisJis-r))dg,. 
Jo 

For the Riemann-Sticltjes integral, we will also give a version of [3], Proposition 2.4. 

Proposition 3.2. Lei g € W^~"'°°(0,T). Assume that a satisfies (HI) and[s,t] C [0,T]. 
/ e Wo'^^i-r, T- M'^), then 

Gi''\f) e Ci-"(s,t;R'^) C W;^-'^{s,t;R'') 

and 

1- l|G^"^(/)lli-„(M)<A„(5)d(3)(l + ||/|U,oo(-.,t-.)), 

nG\^\f)\LMs,t) < + ll/IU.A(-.,*-.)) 

for all A > 1, where d^, i e {3,4}, are positive constants independent of X, f and g. 

Finally, we recall [3], Proposition 2.6. Consider (p{j,a) defined such that (p{j,a) = 2a 
if 7 = f , (p(7, a) > f + if 1^ < 7 < 1 and ifij, a) = a if < 7 < Note that 

(fiij, a) £ [a, 2a]. 

Proposition 3.3. Let g e W^^"'"" {0,T) . Assume that a satisfies (HI) and (H3). // 

/£ Wo"'°°(-r,T;M'*), then 

l|G(-H/)L,.<Ao(ff)rf^^^(l + |^5S) 
for all A > 1, where d^^^ is a positive constant depending only on a, f3,T,d,m and Bq a = 

II&oIIlvc. 

4. Deterministic integral equations 

In this section, we give all the deterministic results. 

For simplicity, let us assume that T = Mr. Set < a < G W^~°''°° {0,T;R'') and 
77 G W(^'°°(— r, 0;R5|_). Consider the deterministic stochastic differential equation on 

x{t)^7iO)+ f b{s,x)ds+ f a{s,xis-r))dg, + y{t), tG(0,T], 
Jo Jo 

(4.1) 

x{t)=r]it), te[-~r,0], 



where, for each i = 1, . . . , d. 



2/'(i)= max(zXs))-, te[0,T], 

s<£[0,t] 
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and 

z{t)=r]{0)+ [ b{s,x)ds+ [ a{s,x{s-r))dgs, te[0,T]. 
Jo Jo 

The existence and uniqueness result reads as follows. 

Theorem 4.1. Assume that b and a satisfy hypotheses (HI) and (H2), respectively, 
with p= 1/a and < a < min{i,/3}. The equation (4-1) then has a unique solution 
Wo"'°°(-r,r;R$). 

Proof. To prove that equation (4.1) admits a unique solution on [— r, T], we shall use 
an induction argument. We will prove that if equation (4.1) admits a unique solution 
on [— r, nr], then we can further prove that there is a unique solution on the interval 
[—r, {n + l)r]. 

Our induction hypothesis, for k < M, is the following: 
(Hk) The equation 

xHt)^v{0)+ f b{s,x^)ds+ f a{s,x''{s-r))dgs+ykit), te[0,kr], 
Jo Jo 

x'=(t) = r/(t), te[-r,0], 
where, for each i = 1, . . . ,d, 

max(4(s))-, te[0,kr], 

sG[0,i] 

with 

Zkit)^r]{0)+ f b{s,x'')ds+ f a{s,x''{s-r))dg„ te[0,kr], 
Jo Jo 

has a unique solution x^ G Wq'°" {~r,kr\W^) . 

The initial case can be easily checked. Assume now that (Hi) is true for all i < n, where 
n < M. We wish to check (Hn+i). 

Clearly, for t G [—r,nr], a;"+^(t) will coincide with x"{t), the solution of the equation 
of (Hn). Moreover, for [—r,nr], yn+i{t) will coincide with i/„(t). We can therefore 
write the equation of (Hn+i) as 



X 



n+li 



{t)^7]{0)+ b{s,x"+^)ds+ cr{s,x''{s-r))dg. 



+ y„+i(t), tG [0,(n + l)r], (4.2) 
x-+\t)=r,it), <G[-r,0]. 
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Moreover, using the notation introduced in the previous section, we have 

it) = 77(0) + F^''^ + G(^) (x") + ijn+i {t), te [0, {n + l)r] , 

x-+\t)^i^{t), te[-r,0]. 

The proof wih be divided into three steps: 

1. if x""*"^ is a solution of (Hn+i) in the space C(— r, (n+ l)r;]]i5|), then x""''^ S 
<'°°(-r,(n + l)r;Rt); 

2. the solution is unique in the space C(-r, (n+l)r;M^); 

3. there exists a solution in the space C(— r, {n + l)r;R'^). 

Step 1: If x^^^ is a solution of (H„-|_i) in the space C(— r, + l)r;M!^), then x^^^ G 
Wo"'°°(-r, (n + l)r;]R^). Wc can write 

ll*^ ||a,oo{— r,(n+l)r) 



te[-r,(n+l)r] V J-r ~ '^j" 

< sup [\x^t)\+ ^i±^^ds 



* |x"+i(0-x"+i(s)| 



te[nr,(n+l)7-] 



— 11*^ II q,cxd( — r.nr) H*^ || Q,cxD(nr, (n+l)r) 

/■"'■ -x"(nr) + x"(rir) -x"(s)| , 
+ sup / — — as 

te[nr,{n+l)r]J-r ~ ■Sj" 

»+in , \x^+\t)~x^^{nr)\ 



— ^Il*^ ||q,oo( — r.nr) ^~ 11-^ II a, c«(nr, (n+l)r) SUp 



tg[„r,(n+i)r] a{t-nr)'^ 

= 2||2;"||„^oo(_r,„r) + ^1 + ^2, 

where 



Ai = ||x"+i| 



Q,oo(nr,(n+l)r) i 

|a;"+i(t) -a:"(nr)| 



A2 = sup 

te[nr,(>i+i)r] a(t - nr)°' 

From our hypothesis, it is clear that oo(_r „r) < oo. So, to complete this step of 

the proof, it suffices to check that < oo and ^2 < oo. 
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We begin with the study of Ai. Clearly, 

Ai<|r;(0)| + ||^^W(z"+i)IU,^(„,,(„^i),) 

(4.4) 

+ l|G'r'^H2;")IU^oo(nr,(n+l)r) + II ^"+1 (") 1 1 a,oo(nr,(n+l)r) • 

One of the keys to our proof is the study of the behavior of y. We note that from its 
definition, it is clear that if is increasing at t (i.e., y'^{t) > y'^{t — e) for e < Eq small 
enough), then 2/'(i) = —z'^{t) > 0. Moreover, by construction, y^{s) > —z^{s) for any s. 
So, if is increasing at t, then for all s <t, 

\y\t) - y\s)\ = y\t) - y\s) < -z^t) + z\s) < \z\t) - z\s)\. 

For t G {nr, [n + l)r) and i G {1, . . . , d}, set 

tl = mi{u-y\u)=y\t)}\Jnr. 

Since is increasing, we note that y*(s) = 2/*(io) fo^ s G [t}),t]. Then, 

* \vl+i{t)-yl+i{s)[^^^ f'° \y'n+i{th)-yn+ii^)\ 



[t-s] 



a+1 



nr 



(t- 


s)«+i 


Wn+lit'o) 


-2/r\+l(5)l 




-s)"+i 


^n+l(^o) 


-4+i(s)l 



< r^^^^±^^?^^-T^^^d. (4.5) 



f*0 

(i^-s)"+i 



On the other hand, we have 



K+i{t)\^\yU^m< sup K+,{s)\. (4.6) 

0<s<tj| 

So, combining (4.5) and (4.6), we have that 

llyn+ill 

< d(||z„+i|| 

a,oo(nr,(n+l)r) ~f~ II -^n || oo(0, nr )), (4.7) 

where we can use the bound 

INn+lll Q:,oo(nr,(n+l)r) 

< |rKO)| + WF^'H^-'-'n^^^inrAn+Dr) + 1 1 (x" ) I L,„o(„.,(„+l).) ■ 

Now combining (4.4), (4.7) and (4.8), we obtain that 

< (d + l)(h(0)| + ||fW(x"+1)|U^^(„,^(„+,),) 

+ \\Gi"\xnL^ooinrXn+l)r))+d\\Zn\UiO,nr)- 



(4.8) 



(4.9) 
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From our hypothesis and Propositions 3.1 and 3.2, it is easy to obtain that H^^n !|oo(o, nr) < 
oo. It therefore only remains to check the norms of the Lebesgue and Riemann-Stieltjes 
integrals. 

On the one hand, 

a,oo(nr,(n+l)r) 

< sup / \bis,x"+')\ds+ ■'■^ \) ds 



tG[nr,(n+l)r] \Jo 



(t^s) 



< sup { [Lo sup |x"+i(M)|+foo(s) ds (4.10) 

Lo X*(sup„^<^<„ \x"+^{v)\ + bo{u)) du 
^ — n as 

,.1-2q 



< Lo ( r+ ^ ) ||.T"+ioo(-.,(„+i).) + T'-" + I|6o||li/o 



On the other hand, to study the Young integral, we will use Proposition 3.2 and the fact 
that x" £ WQ'°°{-r,nr]R^). For any A > 1, 

l|G^'^n^")IL.oo(„..(„+i).) < e'<"+^)n|G(^Hx")L,,(_(„+i),) 

(4.11) 

<^^ti^e^^"^'^'^(i + ||x"|U..(_,„.)). 

So, combining (4.9), (4.10) and (4.11), and using the facts that ||-i;"||Q,A(-r,Tir) < °o and 
ll^;""''"'" l|oo(-r,(n+i)r) < OO for A > 1, wc get that Ai < oo. 
We now deal with the term A2. We can write the decomposition 

\x"+'it)-x"{nr)\ ^ \Jlb{s,x"+')ds\ ^ | /^^ a(s, x"(s - r)) dg,| 
(t-nr)°' ~ (t-nr)"' (t - nrY 

(4.12) 

^ \yn+i{t) - yn+i{nr)\ 
[t - nr)" 

Using the same arguments as in (4.5), we get 
IVn+iit) - yn+i{nr)\ 

t£[nr,{n+l)r] [I — nr ) 

<d sup (4^13) 

te[nr,{n+l)r] [t — nr)'^ 

( \jlb{s^c-+^)ds\ |X*^a(s,a:"(s-r))dgJ\ 
< a sup — h sup 



nr 

[rir,(ri + l)r] [t — Tir)"' te[nr,{n+l)r] {t — UtY 
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Moreover, from Propositions 3.1 and 3.2, we obtain the estimates 



|/„^&(s,a;"+i)ds| \\p(b)( n+i.u i- 

t£lnr,{n+l)r] [t — Tir ) 



2a 



-l2^(l + ll^"^ IU(-r,(n+l)r)), 



(4.14) 



sup Ci-nr^" -11^ \^ )\\l-a{nr,(n+l)rf 

te[nr,(n+l)r] l*- "■' j 

— y,2a — 1 ^ 11"^ ||a,c>o( — r,nr))- 



(4.15) 



< cx), and combin- 



So, using the facts that ||a;"||Q,.oo(-r.nr) < oo and ||a;"^"'" ||oo(-r,(n+i)r) 
ing (4.12), (4.13), (4.14) and (4.15), we obtain that < oo. 
The proof of the first step is now complete. 

Step 2: Uniqueness of the solution in the space Ci-r, {n + l)r;Wl). 
Let X and x' be two solutions of (4.2) in the space C{—r, {n + l)r;R'^) and choose N 
large enough so that ||a;||oo(-r,(ri+i)r) < N and ||a;'||oo(-r,(ri+i)r) < ^• 
For any te [0, (n + l)r], 

sup \x{s) — x' {s)\ < sup \z{s) — z' {s)\ + sup \y{s) — y'{s)\. 
s6[o,t] se[o,t] se[o,t] 

Moreover, using Lemma A. 2, we have 

sup \y{s)-y'{s)\<Ki sup \z{t) - z'{t)\. 
sG[OJ] se[o,t] 

So, combining the last two inequalities, we get that 

sup |x(s) — a;'(s)| < (1 + A';) sup \z{s)-z'{s)\ 
se[o,t] se[o,t] 



< {l + Ki) sup 

se[o,t] 



< {1 + Ki)Ln sup 

se[0,t] 



(6(u, x) — b{u, x')) du 



sup \x{v) — x'{v) \ du 

0<ti<ii 



<Li\[{l + Ki) / sup \x{v) — x' {v)\du. 

Now applying Gronwall's inequality, we have that for all t G [0, (n + l)r], 

sup \x{s) — x'{s)\ = 0. 

sG[0,t] 
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So, 



- a;'||oo(-r,(n+l)r) = 



and the uniqueness has been proven. 

Step 3: Existence of a solution in C{—r, (n+ l)r';K'^). 

In the space C(— r, (n + l)r;]R^[), we can deal with the reflection term using the Sko- 
rokhod mapping. Nevertheless, since the coefficient b is only locally Lipschitz, we will 
need to use a fixed point argument in C(— r, (n + l)r;]R'|_) based on Lemma A.l. 

Let us consider the operator 

C : C(-r, (n + l)r; M^) C(-r, (n + l)r; M|) 

which is such that 

C{u){t)=i]{0)+ f b{s,u)ds+ f a(s,a:"(s-r))d5,+2/„+i,,(t), t G [0, (n + l)r], 
Jo Jo 

C{u){t)=rjit), te[-r,0], 
where x" is the solution obtained in (Hn) and if 

Zn+i,u{t)=v{0)+ [ b{s,u)ds+ [ <T{s,x"{s-r))dgs, 
Jo Jo 

then = maxs6[0,t]«+i,„(s))- for alH = l,...,d. 

Note that C is well deffiied. Moreover, if u = C{u), then u{t) = a;"(i) for any t G (— r, nr). 
We will use the notation u* = C{u). 

We next introduce a new norm in the space C(— r, [n + l)r; W^): for any A > 1, 

ll/l|oo,A(-r.(n+l)r) := SUp e^^^\f{t)\. 

te[-r,(n+l)r] 

These norms are equivalent to ||/||c>o(— 7-,(Ti+i)r)- 

We now check that we can apply Lemma A.l. Note first that 



\U ||cxj,A(-r,(Ti+l)r) 

< lk/l|oo,A(-r,o) + |77(0)| + sup e" 

te[0,(n+l)r] 



-At 



sup e 

tG[0,(n+l)r] 



-At 



(t(s,x"(s - r))dgs 



b{s, u) ds 



d sup e-^*|z„+i,„(0| (4.16) 

te[0,(n+l)r] 



<\\v\\oc,x{-r,o) + id+lMO)\ + {d + l) sup e 

te[0,(ri+l)r] 



-At 



b{s, u) ds 



+ {d+l) sup e~^* 

tG[0,(Ti-fl)r] 



a{s,x"{s ~r))dgs 
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where we have used computations similar to those in (4.6). Indeed, for fixed t, let ti 
mi{u;y'{u)=y'(t)}. Then, 

e-^*|y;+i,„WI <c-^*My;+i,„(ii)l <e-^*M4+i,Jii)l, 
and taking suprema, we have 



sup e ^*\yn+i,u{t)\<d sup e ^*|z„+i_„(<)| 

te[0,{n+l)r] te[0,(n+l)r] 



(4.17) 



Moreover, we have 



sup e 

te[0,(n+l)r] 



-Ai 



6(s, u) ds 



-At 



< Lq sup e' 

te[0,(n+l)r] 



< ^||w||oo,A(-r,(n+l)r) + ^ II ^0 || LVo 



/ sup |M(z;)|ds+ sup e / bo{s)d 

Jo -r<v<s tg[o (n+l)?-! Jo 



(4.18) 



A 

and, from Proposition 3.2, 



sup e 

tG[0,(n+l)r] 



-At 



(T(s,a;"(s-r))dgs 







< ^ff}f\ ^+\\^''\UH-r,nr)). (4.19) 



So, combining (4.16), (4.18) and (4.19), we have 

k*||oo,A(-r,(n-|-l)r) < Mi{X) + A/2 ( A) || W || oo, A(-r,(n-|-l)r) , 

where 

Mi(A) = MU^^r.o)) + id+lMO)\ + ^^^t¥^II&o||lVo 



M2{X) = {d+l)Lo 



Choose A = Ao large enough so that M2(Ao) < ^ . Then, if \\u\\^ Xo{~r,{n+i)r) < 2Afi(Ao), 
we have 

ll"*l|oo,Ao(-r-,(n+l)r) < 2A/i(Ao) 

and so C{Bo) C Bq, where 

B„ = {ue C(-r, {n + l)r; R^); ||u|U,;,„(_,,(„+i),) < 2Afi(Ao)}. 
The first hypothesis in Lemma A.l is thus satisfied with the metric po associated with 

the norm || • ||oo,Ao(-r,(n-f l)r)- 
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To complete the proof, we only need to find a metric pi satisfying the second hypothesis 
of Lemma A.l. 

Note first that if u e Bq, then ||-u||oc,(_r,(n+i)r) < 2e-^o("+i)''Mi(Ao) := iVo- Consider 
u,u' Bq and A > 1. We then have 

\\C{u) - C{u')\\oo,X{-r.{n+l)r) < SUp C"^* |2;„+i,„ (t) - z'^^^j ^„ (t) | 

te[0,(n+l)r] 

+ sup e"^'|?/„+i,„(t) 

te[0,(n+l)r] 

From Lemma A. 2, note that given t G [0, {n + l)r], there exists some t2 <t such that 

\yn+l,u{t) - y'n+l^ui't)\ < Ki\Zn+l,u{t2) {t2)\. 

So, 

e-^'\yn+iAt)-yUiAt)\<Kie-^'-\zn+iAt2)~z'„^,^^{t2)\ 
and it easily follows that 

sup e-^*\yn+i,uit) -yn+i,uit)\<Ki sup e-^*|z„+i,„(i)-z;+i_„(t)|. 

te[0,(n+l)r] te[0,(n+l)r] 



Then, 



\\C{u) - C{u')\\ac,\(-r,{n+l)r) 

<{l + Ki) sup e~^* / |6(s,w) -&(s,M')|ds 

te[0,(n+l)r] Jo 

<LAro(l + A'i) sup e""^* / sup \u{v) - u' {v)\ ds 
te[o,{n+i)r] Jo a<v<s 

<LNoil + Ki) sup [ e"^(*~^)e-^" sup \u{v) - u' {v)\ ds 

te[0,(n+l)r]Jo -r<v<s 



< LNg{l + A';)-||m - u'||oo.A(-r,(n+l)r)- 



So, if we choose A = Ai such that ^"o^^+-^'^ < 1^ then the second hypothesis is satisfied 
for the metric pi associated with the norm || • ||oo,Ai(-r,(n+i)r) s-nd = a\"*"^^'^ • ^ 
We now check that the solution is (1 — a)-H61der continuous. 

Proposition 4.2. Assume that b and a satisfy hypothesis (HI) and (H2), respectively, 
with p= 1/a and <a< min{i,/3}. The solution x of equation (4-1) then belongs to 
Ci-"(0,T;R'') with 

l|a;||l-a(0,T) <d^'Hl+A„(5))(l+||x|U,oo(-r,T)), 

where d*-^-* is a positive constant independent of f and g. 
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Proof. Note that 

l|-^IIl-«(O.T) < \\z\\l-a{0,T) + ll2/||l-a(0,T)- 

For fixed t G [0, T], set = inf{it < t; y^{u) = Then, is increasing in and it is 

easy to check that — y'^{s)\ < \z\t^) ~ for ali s £ (0,^*). For all s S (0,^*), it 

thus holds that 

\yKt)~yKs)\ ^ \z^{u)-z%s)\ 



{t - s)i-" - [u - sy-" 

and it then follows easily that ||2/||i-a(o,T) 5; '^II^IIi-q(o.t)- So, 

l|a;||l-a(0,T) < (d+ l)||2;||i_„(o,T) 

<(d+l)(h(0)| + ||FW(.T)||,_ 
Using Propositions 3.1 and 3.2, we easily complete the proof. □ 

We will now give an upper bound for the norm of the solution. Recall the definition 
of <y9(7,a): 

(2a, 7 = 1, 

2a - 1 1 - 2a 

1 -2a 

a, 0<7< . 

I. 1 — a 

Lemma 4.3. Assume (HI), (H2) and (H3). T/ie unique solution of equation (4-1) then 
satisfies 

||x|U,ooM<4'^(||r/|U,oo(-..o)+A„(g) + l)exp(T(dW+4'^A„(.g)i/(i-^^^^ 

Proof. First, we need to obtain an upper bound for ||a;||a,A(r)- We begin with the esti- 
mates 

\\x\UMr)<MU.x{-rfi)+ sup e-^' ( \x{t)\ + f Ml:i|Mds 

te[o,T] \ J-r It-sj 

<\\vLM-rfi)+ sup e-^* |x(i)|+ / ^ds 
te[o,T] \ J-r (t-s) 

(4.20) 

'\x{t)-x{s)[^^; 



tmrf J-r (t-sY 



Moreover, 

l|a;|U,A(o,T) < lkllQ,A(o,T) + l|y|Ua(o,T) (4.21) 
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and using the same arguments as in (4.5) and (4.6), we have, 

l|y|U,A(0,T)<rfp|U,A(0,T). (4.22) 

We also know that 

II^IUa(o,T) < mi + \\F^'Hmo.Mo,T) + l|Gl"n^)IU,A(o,T)- (4-23) 
So, combining (4.20)-(4.23) and applying Propositions 3.1 and 3.3, we get that 

l|a^lU,A(r) 



\\x\Ux{r) < hlU,A(-r,0) + {d + lMO)\ + {d + 



2a X^~" 



B 



where 



Bi = sup -— / |6(s,x)|ds, 
tefo.Tl "'o 



te[o,T] Jo 



B2 = sup 
te[o,T] 



-At 



a{s,x{s - r))Ags 



^3= sup -—\y{t)\. 
te[o,T] t 

Using the same arguments as in (4.17), we get that 



sup — |y(t)|<d sup —-\z{t)\<d{Bi+B2 
tG[0,T] £ te[o,T] ^ 



(4.24) 



with 



te[OT] J-r (t-s)"+i 
It remains to investigate B. We can decompose this term as follows: 

w |2;(i)-r;(0)| , ^. 1^7(0) -r;(s)| , 

B< sup c"^* / / /i; ds+ sup c-^' / ' \ i\ ^' ds 

1 e--^* 

<- sup -— |x(t) -7y(0)| + ||77|U,A(-r.O) (4.25) 

ate[o,T] ^ 

<i(Bl+B2+S3) + ||77|U,A(-r,0), 
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We now consider Bi and B2 ■ For Bi , we can write 

Bi < sup / ILq sup \x{u)\ + bo{s)j ds 

te[0,T] Jo ^ -r<u<s ' 

_ g-A(t-s) g-At 

<io sup e '*|a;(s)| sup / — ds + sup ||&o||li/o 

< ioA"-ir(l - a)||x|U,;,(,) +C„A2"-1||&o|Li/=- 
Next, we obtain a bound for We will use the hypothesis (H3). 



Bi< sup -— A„(g)( / ds 

te[o,T] I 



- a 



s 

/■* r |cr(s,x(s ~ r)) ~ a{y,x{y - r))| \ 
'„ io ^^^F^^ ^'"V 



^ sup — A„(g) Ao/ ds 



,r / /'^ / |a;(s — r) — xfw — r)| 1 \ , , 



< sup Aa{g)(Ko- e-^* + Xo— / r^4^ds 

tG[o,T] V 7-r [s + r)" 

+ aMo||x|U.My_^ ^^3^ds+^^^— ^j. 

Now, using the inequalities 

sup <''e-^* < fYVe"", 

/ 7 ^ ds = e-^-^ / du 

J-r (^-s)" Jo (t-u + r)" 

<e-y^^^d.<e-A-r(l-a), 
the Holder inequality and the fact that \ f{s)\'' < |/(s)| + 1, we get 

e:^ f^-^ ^(sr < g-A.,.(.,.).-2.+i-. 1^(^)1 ds^ ' 

t" (s + r)" - Vi-r (s + r)"^(T.") y 

< C„ ^ Te-^* f 1 + /*" , '^^'^j , ds) 
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< C'a,7,T( 1 + l|-i;|U,A(r)e / ^— -7— ^ du 



r-t g-A(t-«) 

<C„,^,t(1 + ||x|U,,(,)C-^^A^(^'")-1), 
where we have used the fact that ip{j, a)j — 2a + 1 — 7 > 0. We thus obtain 



^2<A„(g)| 



\l-a\ c J ^ (/3-a)(/?-a + l) 

<A„(g)C„,^,^(l + A2"-i+c-^'-A^(^'")-i||z|U,;,(,)). 
Finally, combining (4.24), (4.25) and the estimates for B,Bi and B2, we have 

||a;|U,A(r) <Mi(A)+M2(A)||a;|U,A(,) 

with 



A/i(A) = 2llrylU,;,(_,,o) + (^ + 1) ( h(0)| + A„(5)d(5) + 

^(d(2) + ||6olLv.+A.(g))Y 



^'^2(A) = §^^(Aa(3)(d(^) + C^,p.,) + LoT{l -a)+ S^^). 

Choosing A = Aq large enough so that ^/2(Ao) = ^, we then have 

IklU.AoM <2A/i(Ao). 

Set 

Ao - [2C„,,(d(2) + A„(g)(d(5) + 1) +io)]'/<'-'^'"^"» 

< d„(2C„,,(d(2) + i^))V(i-^(7,a)) ^ A„(g)V(i-^(7,"))d„(2C„,,(l + d(5)))^/(i-'^(''-")) 
<dW+42)A„(g)i/(i-^(7.")) 

with 

dW=d„(2C„,,(d(2)+Lo))^/(^-^(^'"», 
d(2)=d„(2C„,,(l + d(^)))^/(^-'^('''"». 

This implies that 

l|x|U,ooM <exp(T(dii) +42)A„(5)V(i-^(7.a))))2M,(Ao) 
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and the proof is then easily completed. Note that we can choose da\da'^ which do not 
depend on /3 or 7. □ 

5. Stochastic integral equations 

In this section, we apply the deterministic results in order to prove the main theorem of 
this paper. 

The stochastic integral appearing throughout this paper, u{s)dWs, is a pathwise 
Riemann-Sticltjes integral and it is well know that this integral exists if the process u{s) 
has Holder continuous trajectories of order larger than 1 — H . 

Set a e (1 — i7, i). For any 5 G (0,2), by Fcrnique's theorem (see [1], Theorem 1.3.2), 
we have 

£;(exp(A„(M^)*))<oo. 

Then, if u = {wt,t G [0,T]} is a stochastic process whose trajectories belong to the space 
W^'^ (0, T) , it follows that the Riemann-Sticltjes integral /J" u{s) dWs exists and we have 
that 

i-T 



u{s) dWs 



<G\\u\\ 



where G is a random variable with moments of all orders (see [11], Lemma 7.5). Moreover, 

if the trajectories of u belong to W(^'°°(0,T), then the indefinite integral /^^ u{s)dWs is 
Holder continuous of order \ ~ a and with trajectories in W(5*'°°(G,T). 

The existence and uniqueness of a solution then follows from Theorem 4.1. In order 
to get the existence of a moment of any order, we need only note that if a < (2 — 7)/4, 
then l/(l-(y9(7,a)) <2 and £;(exp(CA„(iy)i/(i-^('''")))) < cx). 



Appendix 

In this appendix, wc just give a fixed point theorem, well posed to our problem, and 
recall a result with some properties of the solution of the Skorokhod problem. 

Lemma A.l. Let {X^p) he a complete metric space, and po and pi two metrics on X 
which are equivalent to p. If C:X ^X satisfies: 

1. there exists some tq > 0, xq G X such that if Bq = {x £ X; po{xQ,x) < tq} then 
CiBo) C Bo; 

2. there exists some a G (0, 1) such that pi{C{x) , C{y)) < api{x,y) for all x,?/ G Bq, 
then there exists some x* G C{Bq) C X such that x* ~ Cix*). 

Proof. For all ?i G N, set 

Xn-\-\ — ^{Xn). 
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Clearly, x„ £ C{Bo) for all n g N. Moreover, 
and 

) H [-pi(x„+i,x„) 

a" 

< a"(aP"^ + ••■ + «+ l)pi(xi,a:o) < pi(xi,xo) -J- 

1 — a 

as n — ^ oo . 

Since {X, p) is a complete metric space and Bq is closed in X, there exists some x* G Bq 
such that *. Furthermore, from the second hypothesis of the lemma, we get that 

Pl{C{Xn),C{x*)) < api{Xn,X*). 

Since pi(x„,x*) -> 0, £(.t„) ^> C{x*) and it follows that x* = C{x*). □ 

Lemma A. 2. For each path z G C(R+,R''), there exists a unique solution {x,y) to the 
Skorokhod problem for z. Thus, there exists a pair of functions (0, 1^9) :C+(M+,M'') — >■ 
C+(IR+,R^'') defined by {(l){z),ip{z)) = {x,y) such that the following holds: There exists 
a constant Ki > such that for any zi, Z2 G C+(R+, K''), we have, for each t>0, 

-(?!)(z2)||oo,[0,i] < A';||zi -Z2||oo,[0,t], 
\\ip{zi) - ip{z2)\\oo,[0,t] < -^i^zll^l -^:2l|oo,[0,t]- 

Proof. See [6], Proposition A.0.1. □ 
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